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Abstract —  This  paper  presents  regression  method’s  in 
estimating  the  free  vibration  analysis  and  compared  with 
SDSST  method.  In  this  study,  the  free  vibration  analysis 
of  the  cross-ply  laminated  composite  cylindrical  shallow 
shells  has  been  studied  using  shear  deformation  shallow 
shell  theory  (SDSST).  First,  the  kinematic  relations  of 
strains  and  deformation  are  given.  Then,  using 
Hamilton ’s  principle,  governing  differential  equations  are 
developed  for  a  general  curved  shell.  Finally,  the  stress- 
strain  relation  for  the  laminated,  cross-ply  composite 
shells  are  obtained.  By  using  some  simplifications  and 
assuming  Fourier  series  as  a  displacement  field,  the 
governed  differential  equations  are  solved  by  the  matrix 
algebra  for  shallow  shells.  Employing  the  computer 
algebra  system  called  MATHEMATICA;  a  computer 
program  has  been  prepared  for  the  solution  [1].  The 
results  obtained  by  this  solution  are  compared  with  the 
results  obtained  by  (ANSYS)  programs.  In  this  article, 
regression  method’s  and  SDSST  method’s  abilities  in 
estimating  the  free  vibration  with  the  laminate  number, 
aspect  ratio,  thickness  ratio,  curvature  ratio  and 
orthotropic  ratio  variables,  are  compared  with  different 
and  similar  aspects.  In  comparing  with  linear, 
interaction,  quadratic  and  pure  quadratic  models,  which 
are  constructed  with  multiple  linear  regression  approach, 
the  quadratic  model  provides  better  results. 

Keywords — Anisotropy,  Finite  Element  Method  (FEM), 
Multi  Linear  Regression,  Shell  Theory,  Structural 
Composites,  Free  Vibration. 

I.  INTRODUCTION 

A  structural  composite  material  consists  of  two  or  more 
constituents  combined  on  a  macroscopic  scale  to  form  a 
useful  material.  Different  materials  must  be  put  together 
in  a  three  dimensional  body.  The  goal  of  this  three 
dimensional  composition  is  to  obtain  a  property  which 
none  of  the  constituents  possesses.  In  other  words,  the 
target  is  to  produce  a  material  that  possesses  higher 
performance  properties  than  its  constituent  parts  for  a 
particular  purpose.  Some  of  these  properties  are 


mechanical  strength,  corrosion  resistance,  high 
temperature  resistance,  heat  conductibility,  stiffness, 
lightness  and  appearance.  In  accordance  with  this 
definition,  the  following  conditions  must  be  satisfied  by 
the  composite  material.  It  must  be  manmade  and 
unnatural.  It  must  comprise  of  at  least  two  different 
materials  with  different  chemical  components  separated 
by  distinct  interfaces.  It  must  possess  properties,  which 
none  of  the  constituents  possesses  alone  and  that  must  be 
the  aim  of  its  production.  The  material  must  behave  as  a 
whole,  i.e.  the  fiber  and  the  matrix  material  (material 
surrounding  the  fibers)  must  be  perfectly  bonded. 
Structures  composed  of  composite  materials  offer  lower 
weight  and  higher  strength  and  stiffness  than  those 
composed  of  most  metallic  materials  [2] . 

Shells  are  common  structural  elements  in  many 
engineering  structures,  including  concrete  roofs,  exteriors 
of  rockets,  ship  hulls,  automobile  tires,  containers  of 
liquids,  oil  tanks,  pipes,  aerospace  etc.  A  shell  can  be 
defined  as  a  curved,  thin- walled  structure.  It  can  be  made 
from  a  single  layer  or  multilayer  of  isotropic  or 
anisotropic  materials.  Shells  can  be  classified  according 
to  their  curvatures.  Shallow  shells  are  defined  as  shells 
that  have  rise  of  not  more  than  one  fifth  of  the  smallest 
planform  dimension  of  the  shell  [2].  Shells  are  three- 
dimensional  (3D)  bodies  bounded  by  two  relatively  close, 
curved  surfaces.  Since  the  3D  equations  of  elasticity  are 
complicated,  all  shell  theories  (thin,  thick,  shallow  and 
deep,  etc.)  reduce  the  3D  elasticity  problem  into  a  2D  one. 
This  is  done  usually  by  Classical  Lamination  Theory- 
CLT  and  Kirchhoff  hypothesis.  A  number  of  theories 
exist  for  layered  shells.  Many  of  these  theories  were 
developed  originally  for  thin  shells  and  based  on  the 
Kirchhoff-Love  kinematic  hypothesis  that  straight  lines 
normal  to  the  undeformed  mid- surface  remain  straight 
and  normal  to  the  middle  surface  after  deformation.  Many 
studies  have  been  performed  on  characteristics  of  shallow 
shells.  Qatu  [2]  uses  energy  functional  to  develop 
equation  of  motion.  Reddy  [3]  are  presented  effect  of 
shear  deformation  for  laminated  composite  plates  and 


www.ijaers.com 


Page  |  114 


International  Journal  of  Advanced  Engineering  Research  and  Science  (IJAERS)  [Vol-4,  Issue-12,  Dec-  2017] 

https:/ /dx.  doi.  org/1 0. 221 61/iiaers. 4.12.19  ISSN:  2349-6495(P)  /  2456-1908(0) 


shells.  Dogan  [4]  studied  the  effects  of  anisotropy  and 
curvature  on  free  vibration  characteristics  of  laminated 
composite  cylindrical  shallow  shells.  Also  Dogan  [5]  are 
presented  the  mode-shape  analysis  of  the  cross-ply 
laminated  composite  cylindrical  shallow  shells  using 
SDSST  and  finite  element  method. 

There  has  been  a  lot  of  number  of  studies  on  the 
regression  analysis  of  composites.  Dong  etc  [6]  studied 
dimension  variation  prediction  for  composites  with  finite 
element  analysis  and  regression  modeling.  The  studied 
provides  information  to  develop  practical  and  forward- 
looking  dimension  control  techniques  for  composite 
material.  Lee  etc.  [7]  were  studied  on  regression  of  the 
response  surface  of  laminated  composite  structures.  In 
this  study,  the  response  surface  of  composite  laminated 
structures  was  calculated  using  regression  analysis.  Ply 
angles  of  composite  structure  was  investigated  using 
regression  analysis  method.  Satapathy  etc.  [8]  studied 
targeted  material  design  of  fly  ash  filled  composites  for 
friction  braking  application  by  non-linear  regression 
technique.  The  study  demonstrated  non-linear  regression 
technique  successfully  analyzed  materials  design  of 
heterogeneous  composites.  Ziari  etc.  [9]  studied 
predicting  rutting  performance  of  carbon  nano  tube 
asphalt  binders  using  regression  models  and  neural 
networks.  Prediction  models  were  obtained  using 
regression  and  artificial  neural  network  of  rutting 
performance  of  carbon  nano  tube  reinforced  asphalt 
binders  and  performances  of  these  models  were 
compared.  Artificial  neural  network  model  performance  is 
better  than  regression  model  results.  Oladipo  etc.  [10] 
studied  optimization  and  modelling  using  non-  linear 
regression  technique  to  enhanced  removal  of  crystal  violet 
by  low  cost  alginate  and  acid  activated  bentonite 
composite  beads.  In  this  work,  operational  parameters 
investigated  and  for  this  used  non-linear  regression  as 
method  of  analysis.  Yadollahi  etc.  [11]  studied 
application  of  adaptive  network-based  fuzzy  inference 
system  (ANFIS)  and  regression  models  to  predict  the 
compressive  strength  of  geopolymer  composites. 
Materials  properties  of  geopolymer  composites  specified 
very  difficult  for  this  reason  geopolymer  are  widely 
complex  materials.  ANFIS,  linear  and  non-linear 
regression  methods  were  studied  to  determine  best 
approach  of  composite  materials. 

II.  THEORY 

A  lamina  is  produced  with  the  isotropic  homogenous 
fibers  and  matrix  material.  Any  point  on  a  fiber,  and/or  on 
matrix  and/or  on  matrix-fiber  interface  has  crucial  effect 
on  the  stiffness  of  the  lamina.  Due  to  the  big  variation  on 
the  properties  of  lamina  from  point  to  point,  macro - 
mechanical  properties  of  lamina  are  determined  based  on 


the  statistical  approach.  According  to  FSDT,  the 
transverse  normal  do  not  cease  perpendicular  to  the  mid¬ 
surface  after  deformation.  It  will  be  assumed  that  the 
deformation  of  the  plates  is  completely  determined  by  the 
displacement  of  its  middle  surface. 

The  theory  of  shallow  shells  can  be  obtained  by  making 
the  following  additional  assumptions  to  thin  and  thick  (or 
shear  deformation)  shell  theories.  It  will  be  assumed  that 
the  deformation  of  the  shells  is  completely  determined  by 
the  displacement  of  its  middle  surface.  The  derivation  of 
equations  of  motion  is  based  on  two  assumptions.  The 
first  assumption  is  that  the  shallow  shell  has  small 
deflections.  The  second  assumption  is  that  the  shallow 
shell  thickness  is  small  compared  to  its  radii  of  curvature. 
Also,  the  radii  of  curvature  are  very  large  compared  to  the 
in-plane  displacement.  Curvature  changes  caused  by  the 
tangential  displacement  component  u  and  v  are  very  small 
in  a  shallow  shell,  in  comparison  with  changes  caused  by 
the  normal  component  w. 


Fig.  1:  Fiber  and  matrix  materials  in  laminated 
composite  shell 


Using  the  given  equation  below  (Eq.l)  nth  layer  lamina 
plate  stress-strain  relationship  can  be  defined  in  lamina 
coordinates, 
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The  displacement  based  on  plate  theory  can  be  written  as 
u(a,  (3,  z)=u0(a,  (3)+z  cpx  (a,  (3) 

v(a,  (3,  z)=vQ(a,  (3)+zcpy(a,  (3)  (2) 
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w(a,  p,  z)=w0(a,  P) 

where  u,  v,  w,  cpx  and  cpy  are  displacements  and  rotations 

in  x,  y,  z  direction,  orderly.  u0,  v0  and  wo  are  mid-plane 
displacements. 
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Equation  of  motion  for  plate  structures  can  be  derived  by 
Hamilton’s  principle 
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in  which  qx,  qy,  qz,  mx,my  are  the  external  forces  and 
moments,  respectively.  U  is  the  strain  energy  defined  as, 
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Solving  equation  5  gives  set  of  equations  called  equations 
of  motion  for  plate  structures.  This  gives  equation  9  in 
simplified  form  as, 
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When  the  shell  has  small  curvature  it  is  referred  to  as  a 
shallow  shell.  Shallow  shells  are  defined  as  shells  that 
have  a  rise  of  not  more  than  l/5th  the  smallest  planform 
dimension  of  the  shell  [2].  It  has  been  widely  accepted 
that  shallow  shell  equations  should  not  be  used  for 
maximum  span  to  minimum  radius  ratio  of  0.5  or  more. 
For  shallow  shells,  Lame  parameters  are  assumed  to  equal 
to  one  (A=B=1). 

The  Navier  type  solution  might  be  implemented  to  thick 
and  thin  plates.  This  type  solution  assumes  that  the 
displacement  section  of  the  plates  can  be  denoted  as  sine 
and  cosine  trigonometric  functions. 
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Assume  a  shell  with  shear  diaphragm  boundaries  on  all 
edges.  For  simply  supported  thick  shells,  boundary 
conditions  can  be  arranged  as  follows: 

Aa=w0=vo  =  Ma=\\f^  =  0  a  =  0,  a 


Np=w  =u0  =  =  0  a  =  0,b 


(12) 
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The  displacement  functions  of  satisfied  the  boundary 
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where  xm=m7i/a,  yn=n7i/b. 

Substituting  the  above  equations  into  the  equation  of 
motion  in  matrix  form, 
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Following  equation  can  be  used  directly  to  find  the 
natural  frequencies  of  free  vibrations.  The  number  of 
terms  that  taken  into  account  in  the  m  and  n  cycle  is  one 
(i.e.  m=l  and  n=l).  Following  equation  can  be  used 
directly  to  find  the  natural  frequencies  of  free  vibrations. 

[K]{A}  +  (a>mny[M]{A}  =  0  (15) 

III.  NUMERICAL  EXAMPLES 

As  an  example,  a  simply  supported  cylindrical  shell 
which  has  a  ratio  of  radius  of  curvature  (ratio  of  shell 
width/shell  radius)  equals  from  0  to  0.1  (0,  0.025,  0.05, 
0.1)  in  one  plane  and  infinite  radius  of  curvature  in  other 
plane,  has  been  considered  (Fig.  3. a).  The  shell,  in  hand, 
has  a  quadrangle  planform  where  the  ratio  of  plan-form 
dimensions  varies  from  1  to  4  (a/b=l,  2,  4).  As  a  material, 
a  laminated  composite  has  been  used  as  [079070°], 
[079079070°],  [0790707079070°],  symmetrical  cross- 


ply  stacking  sequence  (Fig.  3.b).  Ratio  of  modulus  of 
elasticity  (E1/E2)  which  is  the  ratio  of  modulus  of 
elasticity  in  fiber  direction  to  matrix  direction,  has  been 
taken  from  1  to  50  (1,  2,  5,  10,  15,  25,  50).  Effect  of  shell 
thickness  ratio  that  ratio  of  shell  width  to  shell  thickness, 
a/h=100,  50,  20,  10  and  5,  has  been  examined. 


Fig.  2.  a:  Cylindrical  shallow  shell  [12] 


Fig.  2.b:  Layered  sequence  for  cylindrical  shallow  shell 

[12] 


For  each  case,  the  shell  has  been  solved  with  two  theories. 
First  theory  used  in  the  solution  of  composite  laminated 
shallow  shell  is  shear  deformation  shallow  shell  theory 
(SDSST).  The  second  theory  is  the  Finite  element  model 
(FEM).  Entire  structure  is  meshed  by  finite  elements  in 
this  theory.  Then  assuming  a  suitable  displacement  fields 
for  each  meshing  element,  the  behavior  of  the  structure 
has  been  obtained.  In  this  paper,  finite  element  package 
programs,  ANSYS  [12]  have  been  used.  The  structure  is 
meshed  by  25x25  elements  in  ANSYS  model.  A  8-noded 
quadratic  element  is  considered  as  a  meshing  element 
named  as  SHELL99.  The  element  has  100  layers  to  model 
the  composite  materials  used  in  the  structure.  For  each 
layer  geometric  and  material  properties  is  entered  to 
program.  Furthermore,  thicknesses  of  each  layer,  fiber 
orientations  and  stacking  sequence  must  be  entered 
carefully.  Orthotropic  and  lamination  properties  of  the 
problem  could  be  modeled  by  using  this  element. 
Regardless  of  the  point  used,  programs  provided 
geometric  shape  as  seen  in  Fig.  3. 


Fig.  3:  8-point  quadratic  elements  for  ANSYS 
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iv!  MULTIPLE  LINEAR  REGRESSION  METHODS 

FOR  ESTIMATION  OF  FREE  VIBRATION 

The  regression  equations  are  the  formulization  of  the 
relationship  between  dependent  and  independent 
variables.  These  equations  aims  in  this  section  to  estimate 
the  dependent  variable  which  is  the  free  vibration  using 
the  independent  variables  which  are  the  laminate  number, 
aspect  ratio,  thickness  ratio,  curvature  ratio  and 
orthotropic  ratio  variables  at  the  base  of  multiple  linear 
regressions.  Multiple  linear  regression  models  are  built 
from  a  potentially  large  number  of  predictive  terms 
according  to  independent  variables.  Each  model  has 
different  number  of  coefficients.  Coefficients  of  multiple 
linear  regressions  models  are  determined  according  to 
1260  different  Ansys  models  data  for  values  of  free 
vibration  in  this  study.  The  terms  produced  from  the 
variables  are  highly  effective  in  increasing  the  accuracy 
of  estimates.  For  example,  the  number  of  interaction 
terms  increases  exponentially  with  the  number  of 
predictor  variables.  If  there  is  no  theoretical  basis  for 
choosing  the  form  of  a  model,  and  no  assessment  of 
correlations  among  terms,  it  is  possible  to  include 
unnecessary  terms  in  a  model  that  confuse  the 
identification  of  significant  effects.  Multiple  linear 
regression  models  often  take  the  form  of  something  like 
in  (Eq.16) 

y  =  i 30  +  p1*x1  +  p2*x2  +  p3*x3+ p4*x4  +  p5*x5 
+/?6  *  X1  *  x2  +  /?7  *  Xx  *  X3  +  /?8  *  X3  *  x4  +  /?9  *  xx  *  x5 
+/?10  *X2*X3+  *X2*X4  +  p12  *x2*xs 

(16) 

+/?13  *X3*X4  +  /?14  *X3*X5  +  /?15  *X4*X5+  P16  *  Xl 

+Pl7  *  *2  +  Pl8  *  *3  +  Pl9  *X%+P20*xl+£ 

In  this  formula,  xi,  X2,  X3,  X4,  and  X5  represent  independent 
variables,  y  represents  dependent  variable,  pj  (j=l,...,20) 
represent  coefficients  of  regression,  8  represents  error.  A 
response  variable  y  is  modeled  as  a  combination  of 
constant,  linear,  interaction,  and  quadratic  terms  formed 
from  five  predictor  variables  xi,  X2,  X3,  X4,  and  X5. 
Uncontrolled  factors  and  experimental  errors  are  modeled 
by  s.  The  regression  estimates  model  coefficients  (P_j) 
given  on  xi,  X2,  X3,  X4,  X5  and  y. 

With  the  multiple  linear  regression  approach,  linear 
additive  model,  pairwise  interaction  model,  quadratic 
model  and  pure  quadratic  model  are  formed  in  MATLAB 
program  [13]  by  the  data  set.  The  data  set  consisting  of 
1260  results  of  Ansys  models  is  used  to  determine  the 
coefficients  of  the  regressions  formulas.  The  formula  of 
linear  additive  model  based  on  multiple  linear  regression 
is  shown  in  (Eq.17). 

W  =  —35,82584  -  2,66296  *  LN  +  16,98051  *  AR 

(17) 

+0,17772  *  TR  +  30,76336  *  CR  +  0,53356  *  OR 
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In  this  formula,  W,  LN,  AR,  TR,  CR  and  OR  letters  stand 
for  free  vibration,  laminate  number,  aspect  ratio, 
thickness  ratio,  curvature  ratio  and  orthotropic  ratio 
respectively.  Pairwise  interaction  model  based  on 
multiple  linear  regression  has  a  formula  which  is  shown 
in  (Eq.18). 

W  =  20,03485  -  5,31571  *  LN  +  1,51,42  *  AR 
-0,336  *  TR  +  1,39285  *  CR  -  0,81184  *  OR 
+2,12392  *  LN  *  AR  +  0,03255  *  LN  *  TR 
-1,16235  *  LN  *  CR  +  0,12116  *  LN  *  OR  (18) 
+0,09175  *  AR  *  TR  -  11,44856  *  AR  *  CR 
+0,21859  *  AR  *  OR  +  1,15932  *  TR  *  CR 
+0,00699  *  TR  *  OR  +  1,18130  *  CR  *  OR 
In  this  formula,  W,  LN,  AR,  TR,  CR  and  OR  letters  stand 
for  free  vibration,  laminate  number,  aspect  ratio, 
thickness  ratio,  curvature  ratio  and  orthotropic  ratio 
respectively.  Quadratic  model  based  on  multiple  linear 
regression  has  a  lot  of  terms  about  five  independent 
variables.  These  terms  take  part  in  (Eq.19). 

W  =  2051247  -  0,21684  *  LN  -  16,83334  *  AR 
+0,05440  *  TR  -  23,52538  *  CR  -  0,43582  *  OR 
+2,12392  *  LN  *  AR  +  0,03255  *  LN  *  TR 
-1,16235  *  LN  *  CR  +  0,12116  *  LN  *  OR 
+0,09175  *  AR  *  TR  -  11,44856  *  AR  *  CR  (19) 
+0,21859  *  AR  *  OR  +  1,15932  *  TR  *  CR 
+0,00699  *  TR  *  OR  +  1,1813  *  CR  *  OR 
-0,55769  *  LN2  +  3,56684  *  AR2  -  0,0037  *  TR2 
+240,58983  *  CR2  -  0,00746  *  OR2 
In  this  formula,  W,  LN,  AR,  TR,  CR  and  OR  letters  stand 
for  free  vibration,  laminate  number,  aspect  ratio, 

thickness  ratio,  curvature  ratio  and  orthotropic  ratio 
respectively.  The  formula  of  pure  quadratic  model  based 
on  multiple  linear  regression  is  shown  in  (Eq.20). 

W  =  -35,34822  +  7,76183  *  LN  -  1,36325  *  AR 
+0,56813  *  TR  +  5,84513  *  CR  (20) 

+0,90957  *  OR  -  0,55769  *  LN2  +  3,56684  *  AR2 
-0,0037  *  TR2  +  240,58983  *  CR2  -  0,00746  *  OR2 
In  this  formula,  W,  LN,  AR,  TR,  CR  and  OR  letters  stand 
for  free  vibration,  laminate  number,  aspect  ratio, 

thickness  ratio,  curvature  ratio  and  orthotropic  ratio 
respectively. 

After  the  regression  analyses,  the  models  are  compared 
with  the  performance  criteria  which  are  correlation 
coefficient,  mean  percent  error  and  mean  square  error 
used  in  this  study.  These  methods  on  error  comparison  are 
very  useful  for  finding  the  best  models.  The  correlation 
coefficient  is  used  to  understand  the  correlation  between 
the  results  of  the  models  and  the  Ansys  values.  The  mean 
percent  error  criterion  is  calculated  to  take  into 
consideration  the  error  ratio  according  to  each  result  of 
Ansys  which  has  small  or  big  value  of  free  vibration. 
Another  criterion  to  find  the  best  model  in  this  study  is 
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the  mean  square  error  (MSE)  which  was  used  to  evaluate 
the  performances  of  the  models. 

It  is  seen  that  the  quadratic  model  is  the  best  model  to 
estimate  compressive  strength  when  the  multiple  linear 
regression  methods  are  compared  (Table  1).  Performance 
criteria  for  the  estimation  of  the  results  of  quadratic  model 
are  as  follows:  correlation  coefficient  (R)  0.916,  mean 
square  error  62,2995  mean  percent  error  24,23%. 


Table  1:  Comparison  of  Multiple  Linear  Regression 
Methods  for  Estimation  ofW 


R 

MSE 

MPE 

Linear  model 

0,784 

159,1663 

40,65 

Interaction  model 

0,883 

86,4958 

27,54 

Quadratic  model 

0,916 

62,2995 

24,23 

Pure  Quadratic  model 

0,817 

135,4546 

38,10 

V.  COMPARISON  OF  SHEAR  DEFORMATION 

SHALLOW  SHELL  THEORY  AND 
REGRESSION  WITH  MODELS 

When  is  the  SDSST  model  compared  with  the  quadratic 
model  based  on  multiple  linear  regression  approach,  on 
the  whole  data  set.  MSE  of  the  SDSST  is  less  than 
Regression  Model  (Table  2).  Likewise,  while  the 
Regression  model  has  0.916  R  value,  the  SDSST  model 
has  0.989.  When  we  look  at  the  MPE  value,  the  SDSST 
model  again  has  a  better  value  of  3,75%  compared  to 
24,23%  in  the  Regression  model.  Figure  5  shows  the 
results  of  the  Ansys  and  SDSST  models  in  the  graphical 
environment,  while  Figure  6  compares  the  results  of  the 
Ansys  and  Regression  results.  While  the  results  of  the 
SDSST  model  only  show  large  deviations  in  moderate 
values,  the  regression  results  show  large  deviations  in 
almost  all  values. 


Table  2:  Comparison  ofR,  MSE  and  MPE  of  SDSST  and 
Regression  Model  Results 


R 

MSE 

MPE  (%) 

Regression 

SDSST 

0,91632 

0,98845 

62,2995 

9,2358 

24,23 

3,75 
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Fig.  4:  Comparison  of  Ansys  and  SDSST  Model  Results 
for  each  of  all  data  set 


200 


Number  of  solition 

Fig.  5:  Comparison  of  Ansys  and  MLR  Model  Results  for 
each  of  all  data  set 

VI.  RESULTS  OF  THE  MODELS 

In  this  article,  regression  methods  and  SDSST  method’s 
abilities  in  estimating  the  free  vibration  with  the  laminate 
number,  aspect  ratio,  thickness  ratio,  curvature  ratio  and 
orthotropic  ratio  variables,  are  compared  with  different 
and  similar  aspects.  In  comparing  with  linear,  interaction, 
quadratic  and  pure  quadratic  models,  which  are 
constructed  with  multiple  linear  regression  approach,  the 
quadratic  model  provides  better  results.  With  developing 
quadratic  regression  model,  MPE  declined  to  24,23%,  R 
rose  to  0,916.  The  MSE  value,  which  is  a  very  important 
criterion  for  comparing  model  results  became  less  than 
62,2995  and  positive  developments  are  seen  in  the  whole 
evaluation  criteria.  Similarly,  from  the  models  that  are 
based  on  the  SDSST  theory  provides  the  best  result.  With 
the  SDSST  model,  the  R  rose  to  0,989,  the  MPE  declined 
to  3,7%,  while  the  important  criterion  of  MSE,  which 
shows  high  performance  in  model  comparisons  declined 
to  9,24.  When  comparing  the  quadratic  model  which  is 
the  best  among  the  regression  based  models  and  the 
SDSST,  the  results  of  the  SDSST  based  model  are 
superior.  In  estimating  the  free  vibration,  SDSST  methods 
show  better  performance  than  the  multi  linear  regression 
methods. 

VII.  CONCLUSION 

According  to  outcomes,  the  under  mentioned  conclusions 
can  be  drawn: 

The  quadratic  model  based  on  multiple  linear  regression 
contains  laminate  number,  aspect  ratio,  thickness  ratio, 
curvature  ratio  and  orthotropic  ratio  variables  and  the 
quadratic  model  was  fitted  all  of  the  data  set.  The 
quadratic  model  has  0,916,  62,2995,  24,23%  R,  MSE  and 
MPE  values,  respectively. 

The  results  of  improved  SDSST  model  compare  with 
quadratic  model,  it  is  seen  that  the  SDSST  model  has 
better  results  than  the  quadratic  model  with  0.989  R  ,  9,24 
MSE  and  3,75  %  MPE  values. 

The  fluctuations  of  data  set  of  the  free  vibration  were  very 
well  reflected  using  SDSST  models  constituted  laminate 
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number,  aspect  ratio,  thickness  ratio,  curvature  ratio  and 
orthotropic  ratio  variables.  Furthermore,  the  SDSST 
models  gave  better  reflection  than  the  multiple  linear 
regression  models. 


[10]Oladipo  A.  A.  and  Gazi  M.  (2014).  Enhanced 
removal  of  crystal  violet  by  low  cost  alginate/acid 
activated  bentonite  composite  beads:  Optimization 
and  modelling  using  non-linear  regression  technique. 
Journal  of  Water  Process  Engineering.  2,  43-52. 


It  has  an  incentive  effect  for  future  studies  to  know  that 
both  of  the  methods,  multiple  linear  regression  with 
quadratic  terms  and  SDSST,  produce  better  results  to 
estimate  free  vibration  using  laminate  number,  aspect 
ratio,  thickness  ratio,  curvature  ratio  and  orthotropic  ratio 
variables. 


[11]  Yadollahi  M.  M.,  Benli  A.  and  Demirboga  R.  (2017). 
Application  of  adaptive  neuro-fuzzy  technique  and 
regression  models  to  predict  the  compressive  strength 
of  geopolymer  composites.  Neural  Comput  &  Applic. 
28,  1453-1461. 


The  results  show  that  the  MLR  method  cannot  achieve 
better  results  than  the  SDSST  method.  To  achieve  better 
results,  different  statistical  methods  or  artificial 
intelligence  techniques  should  be  used.  Models  created 
using  artificial  intelligence  techniques  in  particular  can 
achieve  better  results  than  the  SDSST  model. 


[12]  ANSYS  Inc,  User  manuel  Version:  5.3.  Theory 
reference  manuel  and  ANSYS  element  reference. 
http://www.ansys.com 


[13]MATLAB,  The  MathWorks  Inc.,  Natick,  MA. 
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